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The phenomenom of emerging regular spectral features from random interactions is 
addressed in the context of the vibron model. A mean-field analysis links difi'erent regions 
C , of the parameter space with definite geometric shapes. The results that are, to a large 

extent, obtained in closed analytic form, provide a clear and transparent interpretation 
of the high degree of order that has been observed in numerical studies. 
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Random matrix theory was developed to describe statistical properties of nuclear spectra, such as 
average distributions and fluctuations of peaks in neutron-capture experiments |^,|^. In this approach, 
I the Hamiltonian matrix elements are chosen at random, while keeping some global symmetries, e.g. the 

. matrix should be hermitean, and be invariant under time-reversal, rotations and reflections. Specific 

examples include the Gaussian Orthogonal Ensemble (GOE) of real-symmetric random Hamiltonian 
matrices in which the many-body interactions are uncorrelated, and the Two-Body Random Ensemble 
(TBRE) in which the two-body interactions are taken from a distribution of random numbers ||^ . For two 
o : particles, the two ensembles are identical but for more than two particles, unlike the case of GOE, in TBRE 

2 ■ the many-body matrix elements are correlated. As a consequence, also the energy eigenvalues of states 

with different quantum numbers are strongly correlated, since they arise from the same Hamiltonian. 

The latter aspect was investigated recently in shell model calculations for even-even nuclei in the 
sd shell and the pf shell An analysis of the energy spectra of an ensemble of random two-body 
Hamiltonians showed a remarkable statistical preference for ground states with angular momentum and 
' parity = 0"*", despite the random nature of the two-body matrix elements, both in sign and relative 

magnitude. A similar preponderance of 0"*" ground states was found in an analysis of the Interacting 
Boson Model (IBM) with random interactions In addition, in the IBM evidence was found for 
both vibrational and rotational band structures. According to the conventional ideas in the field, the 
occurrence of L = ground states and the existence of vibrational and rotational bands are due to the 
very specific form of the interactions. Therefore, these unexpected and surprising results have sparked a 
large number of investigations to explain and further explore the properties of random nuclei fl^ . 

The basic ingredients of the numerical simulations, both for the nuclear shell model and for the IBM, 
are the structure of the model space, the ensemble of random Hamiltonians, the order of the interactions 
(one- and two-body), and the global symmetries, i.e. time-reversal, hermiticity and rotation and reflection 
symmetry. The latter three symmetries of the Hamiltonian cannot be modified, since we are studying 
many-body systems whose eigenstates have real energies and good angular momentum and parity. It has 
been shown that the observed spectral order is a robust property that does not depend on the specific 
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choice of the (two-body) ensemble of random interactions |pp[j7|,p^ , the time-reversal symmetry |6j , or 
the restriction of the Hamiltonian to one- and two-body interactions § . These results suggest that that 
an explanation of the origin of the observed regular features has to be sought in the many-body dynamics 
of the model space and/or the general statistical properties of random interactions. 

The purpose of this article is to investigate the origin of the regular features that emerge from random 
interactions in a simple exactly solvable model. Hereto we use the vibron model, which is mathematically 
simpler than the IBM, but exhibits many of the same qualitative features. A preliminary account of 
this work has been published in pGf . In Section 2 we present a review of the vibron model by studying 
a schematic Hamiltonian with an arbitrary strength parameter. In Section 3 we discuss the spectral 
properties of an ensemble of random one- and two-body interactions, which are interpreted in Section 4 
in a mean-field analysis. Finally, in Section 5 we present our summary and conclusions. 



II. THE VIBRON MODEL 

The vibron model was introduced in 1981 to describe the rotational and vibrational excitations of 
diatomic molecules and has also found applications in nuclear cluster models Q and meson spec- 
troscopy [ p^ . In general terms, it provides an algebraic treatment to describe the relative motion in 
two-body problems. The algebraic approach consists in quantizing the relative coordinates and momenta 
with vector boson operators with angular momentum and parity = 1~ 



V2' dr, 



and adding an auxiliary scalar boson with = 0+. The building blocks of the vibron model are then 
given by 

s\pi = -1,0,1). (2) 

Sometimes, the scalar and vector bosons are also called vibrons. The 16 bilinear products of creation and 
annihilation operators are the generators of the Lie algebra of U (4) 

s^s, s^'pf,, p^s, pIp^ (a*, = -1,0,1). (3) 

The Hamiltonian and all other physical operators of interest are expressed in terms of the generators. 
Therefore, the total number of vibrons 

7V = sts-f ^ptp^, (4) 

is a conserved quantity. The presence of the scalar boson makes it possible to consider, in addition to the 
three-dimensional harmonic oscillator, also situations in which the oscillator shells are mixed. In addition 
to the total number of bosons N, the eigenfunctions have good angular momentum L and parity P. For 
a more detailed discussion of the vibron model see lESl and references therein. 



A. A schematic Hamiltonian 

For the study of random interactions, it is convenient to first consider a schematic Hamiltonian which 
contains the basic features of the model [E4| 
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H = -cosxp'^ •P+ {sh^ -p^) {ss~p-p), (5) 

with s = s and = (— The dots denote a scalar product with respect to the rotation group. 
The range of the angle x is that of a full period —tt/2 < x < 37r/2, such that all possible combinations 
of attractive and repulsive interactions are included. 

For X = 0, TT the Hamiltonian has a U{3) dynamic symmetry. The spectrum is that of a three- 
dimensional harmonic oscillator, i.e. a series of oscillators shells with 

n = 0,l,...,7V, 

L = n, n - 2, . . . , 1 or 0. (6) 

The parity of the states is P = (—1)" = (—1)^. The energy spectrum is given by 

E = ±n, (7) 

where the + (— ) sign corresponds to x = (tt). 

For X 7''/2, 37r/2 the Hamiltonian has an 50(4) dynamic symmetry. In this case, the harmonic 
oscillator shells are mixed by the Hamiltonian. The spectrum is that of a deformed or Morse oscillator, 
which consists of a series of rotational bands labeled by 

cr = iV, - 2, . . . , 1 or 0, 

i = 0,l,...,a. (8) 
The corresponding energy spectrum is given by 

= ±^vil-^), (9) 

where the + (— ) sign corresponds to x = 7r/2 (37r/2). The ground state band has v — {N — ct)/2 = for 
X = n/2 and v = [N]/2 for x = 37r/2. 



B. Geometric shapes 

The schematic Hamiltonian of Eq. (||) exhibits various geometric shapes (as well as the phase transitions 
inbetween them) which are relevant for the subsequent studies with random interactions. The connection 
between the vibron model, potential energy surfaces, geometric shapes and phase transitions can be 
investigated by means of standard Hartree-Bose mean-field methods [|4|-g6|. For the vibron model, it is 
convenient to introduce a coherent, or intrinsic, state expressed as a condensate of deformed bosons with 
axial symmetry 

\N,a) ^ (cosas^ + sinapl) |0), (10) 



with < a < tt/2. The potential energy surface is then given by the expectation value of the Hamiltonian 
in the coherent state 

l.E{a)^j^{N,a\H\N,a) 

= cosx sin^ a -I- i sinx cos^ 2a. (11) 
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The equilibrium configuration is characterized by the value a — ao for which the potential energy surface 
attains its minimum value 



The solutions can be divided into three different classes or phases 

ao = -7r/2 < x < •^/4 

cos2ao = cotx 7r/4 < x < 37r/4 (13) 

ao = 7r/2 37r/4 < x < 37r/2 

which correspond to an s-boson or spherical condensate, a deformed condensate, and a p-boson conden- 
sate, respectively (see Fig. The nature of the phase transitions at the critical points Xc = 7r/4, 37r/4 
and 37r/2 can be investigated by examining the Hartree-Bose ground state energy and its derivates. The 
ground state energy itself is a continuous function of x 

r isinx -7r/2<x<7r/4 
-£;(X) = <^ icosx-ff^ 7r/4<x<3V4 (14) 

Lcosx + |sinx 37r/4 < X < 37r/2 

The first derivative of E{x) shows a discontinuity at Xc = 37r/2, and hence the phase transition between 
the spherical, or s-boson, condensate and the p-boson condensate is of first order. The phase transitions 
involving the deformed condensate are of second order, since the second derivative of the ground state 
energy is discontinuous at Xc = 7r/4 and 37r/4. 



C. Rotations 



In the previous section, we investigated the equilibrium configurations of the schematic Hamiltonian 
of Eq. (0). Each one of them corresponds to an intrinsic ground state band \N,ao), whose angular 
momentum content depends on the value of ao. The rotational energies can be obtained by applying 
standard many-body techniques [ psf . 

In the coherent, or intrinsic, state of Eq. (0), the rotational symmetry is spontaneously broken. In 
Random Phase Approximation the corresponding spurious excitations are decoupled from the physical 
ones and lie at zero energy. The collective or rotational energies are then determined by the inertial 
parameter associated with the spurious motion 

i^rot = ^L(L + 1), (15) 

where the moment of inertia X is obtained from the Thouless-Valatin formula. The general procedure is 
described in |2^ and applied to systems of interacting bosons in j26j . For the Hamiltonian of Eq. (||) , we 
find that the moment of inertia is given by 

J= . (16) 

cos 2ao sin x 

The ordering of the rotational energy levels is then determined by the sign of the moment of inertia. In 
the following, we examine each one of the equilibrium configurations in more detail. 

(i) For ao = 0, the equilibrium configuration has spherical symmetry, and hence can only have L = 0. 
The moment of inertia is X = 0. 

(ii) For < ao < tt/2, the equilibrium shape is deformed. The intrinsic state is a condensate of N 
deformed bosons, which are superpositions of monopole and dipole bosons with cos 2ao — cot x and 
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'''/4 < X ^ 37r/4. The ordering of the rotational energy levels L — 0,1, 
of the moment of inertia 

^ ^ A^(sinx-cosx) 
sin X cos X 



, N is determined by the sign 



(17) 



For 7r/4<x^7r/2 the moment of inertia is positive and hence the ground state has angular momentum 
L = 0, whereas for for 7r/2 < x < 37r/4 it is negative corresponding to a ground state with L — N. 

(iii) For ap = 7r/2, the coherent state is a condensate of N dipole or p-bosons. The angular momentum 
content is that of a three-dimensional harmonic oscillator shell with N quanta: L = N, N — 2, . . . ,1 orO 
for N odd or even, respectively. The moment of inertia is 



I =-- 



2N 
sinx' 



(18) 



This equilibrium shape arises for 3tt/A < x ^ 3tt/2. For 37r/4 < x ^ the moment of inertia is negative 
and the ground state has angular momentum L = N. For tt < x ^ 37r/2, the moment of inertia is positive 
and the angular momentum of the ground state is L = for iV even and L = 1 for N odd. 

In summary, the schematic Hamiltonian of Eq. (^ gives rise to three different equilibrium configurations 
or geometric shapes, which are separated by phase transitions. An analysis of the angular momentum 
content of the corresponding condensate combined with the sign of the moment of inertia, yields trans- 
parent results for the ground state angular momentum. The results of Table ^ were obtained by assuming 
a constant probability distribution for x on the interval — 7r/2 < x ^ 37r/2. The ground state is most 
likely to have angular momentum L = 0: in 75% of the cases for N even and in 50% for N odd. In 25% 
of the cases, the ground state has the maximum value of the angular momentum L = N. The only other 
value of the ground state angular momentum is L = 1 in 25% of the cases for N odd. The fluctuation in 
the L ~ Q and L = 1 percentages is due to the contribution of the p-boson condensate. The sum of the 
L — and L ~ 1 percentages is constant (75%) and does not depend on the total number of vibrons N. 

An exact analysis in which the Hamiltonian of Eq. (|^) is diagonalized numerically for different values 
of X) yields identical results for the distribution of the ground state angular momenta as are obtained 
from the mean- field analysis. 



III. RANDOM INTERACTIONS 



In this section, we discuss the properties of the vibron model with random interactions, or more 
precisely, with one- and two-body interactions with random strengths. We consider the most general 
one- and two-body Hamiltonian of the form 



N 



1 



N - 1 



-Ho 



where Hi contains the boson energies 

Hi = es s'^ ■ s - tpp^ ■ p, 
and H2 consists of all possible two-body interactions 

H2 = uo i (st X st)(0) . (S X S)(") + ui (st X p^)(i) • {p X 5)(i) 

+ E CAj(p^Xpt)(A).(pxp)W 



(19) 



(20) 



A=0,2 
1 



+Vq 



2V2 



(ptxpt)(o).(5xS)(o)-|-iJ., 



(21) 
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We have scaled Hi by N and H2 by N{N — 1) in order to remove the N dependence of the matrix 
elements. The seven parameters of the Hamiltonian, altogether denoted by 

(x) = {es,ep,uo,ui,CQ,C2,vo), (22) 

are taken as independent random numbers on a Gaussian distribution 

P(x,) = e-^'/2'^V^^^/2^: (23) 

with zero mean and width a. In this way, the interaction terms are arbitrary and equally likely to be 
attractive or repulsive. The spectral properties of each Hamiltonian are analyzed by exact numerical 
diagonalization. The results discussed in this section are based on 100,000 runs. 

In Fig. H we show the percentages of L = 0, L — 1 and L = N ground states as a function of the total 
number of vibrons N. The vibron model shows a dominance of L = ground states. For even values of 
N the ground state has L = in ~ 71% of the cases, and for odd values in ~ 54% of the cases. Similarly, 
the percentage of ground states with L = 1 shows an oscillation between ~ 1% for even values of N and 
^ 18% for odd values. In ~ 24% of the cases the ground state has the maximum value of the angular 
momentum L — N. 

For the cases with a L = ground state, it is of interest to study the probability distribution of the 
ratio of excitation energies 

(24) 

which constitutes a measure of the spectral properties of the vibron Hamiltonian. The energy ratio R 
has characteristic values of i? = 2 for the vibrational or C/(3) limit (harmonic oscillator, see Eq. @)), and 
i? = 3 in the rotational or 50(4) limit (Morse oscillator, see Eq. (p^)). Figs. ^ and ^ show that, both 
for odd and even values of N, the probability distribution P{R) has two pronounced peaks, one at the 
vibrational value of i? = 2 and one at the rotational value of i? = 3. Moreover, for even values of N there 
is a maximum at i? = 0, which is absent for odd values. 

These numerical results are very similar to those found for the IBM in nuclear physics although 
there are some differences as well. Despite the random nature of the interactions both in sign and 
relative magnitude, the spectral properties show a surprisingly large degree of order. In recent studies, 
the tridiagonal form of the Hamiltonian matrix in the U (3) basis of the vibron model was used to establish 
a connection with random polynomials [p^ . However, in general the Hamiltonian matrix is not of this 
form, and one has to look for alternative methods to understand the origin of these regular properties 
in an analytic and more intuitive way. In the next section, we apply the same mean-field techniques 
that were used in Section 2, to the general one- and two-body vibron Hamiltonian of Eqs. (pj|)-([2l]) with 
random interaction strengths. 

IV. MEAN-FIELD ANALYSIS 

The potential energy surface associated with the general one- and two-body vibron Hamiltonian of 
Eqs. (p^)-(]2l|) is given by its expectation value in the coherent state of Eq. ( p^ ) 

E{a) — 04 sin** + 02 sin^ a + oq. (25) 

The coefficients are linear combinations of the parameters of the Hamiltonian 

^ ^ 1 11 1 

ai=r-x = -uo + ui + -cq + -C2 + —7=vo, 

^ ^ 1 

s - X = -eg + Ep - Uo - ui -=vo, 

V6 

ao ^t- X = es + -uq. (26) 



6 



For random interaction strengths, the trial wave lunction ol Eq. ( p^ ) and the energy surface of Eq. ( p5| ) 
provide information on the distribution of shapes that the model can acquire. The value of ao that 
characterizes the equilibrium configuration of the potential energy surface only depends on the coefficients 
04 and 02. Just as for the schematic Hamiltonian of Eq. the parameter space can be divided into 
different areas according to the three possible equilibrium configurations 



ao = 



sin^ ao 



ao 



—02/204 
tt/2 



02 > 0, 04 > — o2 

— 2o4 < 02 < 
02 < 0, 2o4 < — 02 
04 < — 02 < 



(27) 



In Fig. 1^, the three regions in the 02O4 plane are labeled by I for the s-boson or spherical condensate 
(ao = 0), by II for the deformed condensate (0 < ao < 7r/2), and by III for the p-boson condensate 
(ao = 7''/2). They are separated by the separatrices 02 = 0, 04 > for I-II, 02 — — 204, 04 > for II-III, 
and 02 = — 04, 04 < for III-I. The dashed curve corresponds to the schematic Hamiltonian of Eq. (^), 
and is characterized by 04 ~ sinx and 02 = cosx — sinx with — 7r/2 < x ^ 37r/2. In the previous section, 
we showed that this Hamiltonian exhibits three phase transitions: second order transitions at Xc = 7r/4 
and Xc = 37r/4, and a first order transition at Xc = 37r/2. The intersections of the dashed curve and the 
separatrices occurs exactly at the critical points Xc = 7i'/4, 37r/4 and 37r/2. To study the nature of the 
phase transitions for the case of the general Hamiltonian of Eqs. ( p"^^2l| ) we take an arbitrary ellipse in 
the 0204 plane that encloses the origin as its center. It is straightforward to show that the order of the 
phase transitions does not depend on the orientation nor the eccentricity of the ellipse. In other words, 
the phase transitions are independent of the angle under which the separatrices are crossed. 

The distribution of shapes for an ensemble of Hamiltonians depends on the joint probability distribution 
of the coefficients 04 and 02 which, for the Gaussian distribution P{xi) of Eq. (p3|), is given by a bivariate 
normal distribution 



P{a4, 02) — J Y\ P{^i) &{aA — r ■ x) 5{a2 — s ■ x) 



1 



27rVdet M 



exp 



with 



-i(o4 «.)M-(2 



(28) 



M = 





r 


r 




1 f 28 


-30 \ 




s 


s 




" 18 I -30 


75 I 



(29) 



The vectors f and s are defined in Eq. (|2^) . The probability that the equilibrium shape of an ensemble of 
Hamiltonians is spherical can be obtained by integrating P(o4, 02) over the appropriate range I (02 > 0, 
04 > — 02) 



Pi = J da4da2P{a4,a2) 




= 0.396. 



(30) 



Similarly, the probability for the occurrence of a deformed shape can be derived by integrating P(a4, 02) 
over the area II (— 204 < 02 < 0) 
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1 

= — arctan 
27r 




(31) 



Finally, the probability for finding the third solution, a p-boson condensate, is given by 



P3 



I-P1-P2 



0.388. 



(32) 



The angular momentum of the ground state for each of the equilibrium configurations can be estimated 
by evaluating the moment of inertia. Just as in Section 2, we adopt the Thouless-Valatin prescription, 
which leads to the formula 



The moment of inertia depends in a complicated way on the parameters in the Hamiltonian, both explicitly 
as seen in the denominator of Eq. (^3|) and implicitly through ao, which determines the equilibrium 
configuration. For the schematic Hamiltonian of Eq. (0) , it was possible to obtain a closed expression for 
the moment of inertia, since in this case all properties depend on a single parameter x- However, this 
is not the case for the general one- and two-body Hamiltonian that we are considering here. Instead we 
have to solve the problem numerically. 

In Table || we show the probability distribution of the ground state angular momentum as obtained 
in the mean- field analysis. The results are qualitatively the same as those of Table | for the schematic 
Hamiltonian. There is a statistical preference for L = ground states. This is largely due to the 
occurrence of a spherical shape (whose angular momentum content is just L = 0) for almost 40% of the 
cases (see Eq. (|o|)). The deformed shape yields ground states either with L — for positive values of 
the moment of inertia X > 0, or with L = N ior I < 0. The third solution, the p-boson condensate, gives 
rise to ground states with L = N and, depending whether the total number of vibrons is even or odd, 
to L = or L = 1, respectively. The sum of the L = and L = 1 percentages is constant. In Fig. ^ we 
show the percentages of L = 0, L = 1 and L = N ground states, as a function of the total number of 
vibrons N. As is clear from the results presented in Table |l^, the fluctuations in the percentages of L = 
and L = 1 ground states with N are due to the contribution from the p-boson condensate solution. A 
comparison with Fig. |^ shows that the mean-field results are in excellent agreement with the exact ones. 
The difference observed for the L = N percentage arises from the fact that in the exact calculations for 
approximately 5% of the cases, the value of the ground state angular momentum is different from i = 0, 



In Figs. |7| and ^ we show the contribution of each one of the equilibrium configurations to the probability 
distribution P{R) of the energy ratio R of Eq. ( p^ for = 19 and N = 20, respectively. For both cases, 
the spherical shape (solid line) contributes almost exclusively to the peak at R = 2, and similarly the 
deformed shape (dashed line) to the peak at i? = 3, which confirms the vibrational and rotational 
character of these maxima. The occurrence of a peak at small values of R for A^ = 20 corresponds to 
a level sequence L = 0, 2, 1. It is related to the p-boson condensate solution (dotted line), which has 
angular momenta L = N, N — 2, . . . ,0. The first excited L = 1 state belongs to a different band and has 
a higher excitation energy. For odd values of A^ the p-boson condensate has no state with L = 0, and 
hence the peak at i? = is absent. 



We have investigated the origin of the regular features that have been observed in numerical studies 
of nuclear structure models with random interactions. The observed spectral order is a robust property 



1 = 



2N sin^ ao 



(33) 



4:{N-1) ^t;ocos2ao - i(co - C2)sin2Q!o 



1, A^. 



V. SUMMARY AND CONCLUSIONS 
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that arises from a much larger class of Hamiltonians than is usually thought. It cannot be explained by 
the time-reversal symmetry of the Hamiltonian, the choice of a specific ensemble of random interactions, 
or the restriction to at most two-body interactions. 

In this paper, we have carried out an analysis of the vibron model, which is an exactly solvable model 
to describe the relative motion in two-body problems. A numerical study of the vibron model with 
random interactions shows the emergence of ordered features, such as the dominance of ground states 
with L = and the occurrence of vibrational and rotational band structures. In a mean-field analysis, 
it was found that different regions of the parameter space can be associated with particular intrinsic 
vibrational states, which in turn correspond to definite geometric shapes: a spherical shape (~ 40%), 
a deformed one 20%) and a condensate of dipole bosons (~ 40%). Since the spherical shape only 
has L = 0, and the deformed shape and the p-boson condensate with an even number of bosons N in 
about half the number of cases, one finds an i = ground state in approximately 70% of the cases for N 
even and 50% for N odd. The spherical shape gives rise to the occurrence of vibrational structure, and 
the deformed shape to rotational bands. Qualitatively, these results are very similar to those obtained 
in closed analytic form for a schematic vibron Hamiltonian which interpolates between the harmonic 
oscillator (or U{3) limit) and the Morse oscillator (or 50(4) limit). 

In summary, the present results show that a mean-field analysis provides a clear and transparent 
interpretation of the regular features that have been obtained in numerical studies of the vibron model 
with random interactions. In | [T^ we have applied similar methods to the IBM. Since the structure of the 
model space of the IBM is more complicated than that of the vibron model, the analysis becomes more 
difficult, but the final results are qualitatively the same. The fact that these properties are shared by 
different models, seems to exclude an explanation solely in terms of the angular momentum algebra, the 
connectivity of the model space, or the many-body dynamics of the model, as has been suggested before. 
The present analysis points, at least for systems of interacting bosons, to a more general phenomenon that 
does not depend so much on the details of the angular momentum coupling, but rather on the occurrence 
of definite, robust geometric phases such as spherical and deformed shapes. These shapes are a reflection 
of an intrinsic geometry (or topology) associated to the many-body dynamics of the model space which 
is sampled by the statistical nature of the random interactions, but which is quite independent of them. 

For the nuclear shell model the situation is less clear. Although a large number of investigations 
to explain and further explore the properties of random nuclei have shed light on various aspects of 
the original problem, i.e. the dominance of 0^ ground states, in our opinion, no definite answer is yet 
available, and the full implications for nuclear structure physics are still to be clarified. 
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TABLE I. Probabilities of ground states with L = 0, 1 and A'^, obtained in a mean-field analysis of the vibron 
Hamiltonian of Eq. . 



Shape 




L = 


L = 1 


L = N 




Qo = 


3/8 


3/8 










< Qo < 7r/2 


1/4 


1/8 





1/8 




ao = 7r/2 


3/8 


1/4 





1/8 


N ^2k 




3/8 





1/4 


1/8 


iV = 2fc + 1 


Total 


1 


3/4 





1/4 


N ^2k 




1 


1/2 


1/4 


1/4 


N ^2k + l 


TABLE II. PercentaE 


;es of ground states with L = 0, 


1 and A^, obtained in 


a mean-field analysis 


of the vibron 


Hamiltonian of Eqs. ([L9| 


)-(EiD- 










Shape 




L = 


L = 1 


L = N 




ao = 


39.6% 


39.6% 


0.0% 


0.0% 




< ao < 7r/2 


2L6% 


13.3% 


0.0% 


8.3% 




ao = 71/2 


38.8% 


17.9% 


0.0% 


20.9% 


N ^2k 






0.0% 


17.9% 


20.9% 


N = 2k+l 


Total 


100.0% 


70.8% 


0.0% 


29.2% 


N = 2k 






52.9% 


17.9% 


29.2% 


N = 2k+l 
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FIG. 1. Equilibrium configurations of tfie sdiematic Hamiltonian of Eq. (g) as a function of X- 
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FIG. 2. Percentage of ground states with angular momentum L = (solid line), L = 1 (dashed line) and 
L = N (dotted line) in the vibron model with random one- and two-body interactions obtained for 10 < A'^ < 20 
and 100000 runs. 
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FIG. 3. Probability distribution P{R) of the energy ratio R of Eq. in the vibron model with random one- 
and two-body interactions obtained for A'^ = 19 and 100000 runs. 
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15 




FIG. 5. Equilibrium configurations in the 0204 plane: I) spherical shape (a2 > 0, 04 > —02), II) deformed 
shape ( — 2a4 < a2 < 0) and III) p-boson condensate. The dashed curve corresponds to the schematic Hamiltonian 
of Eq. and is characterized by 04 = sinx and 02 = cosx — sinx with 7r/2 < x < 37r/2. 
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FIG. 6. Percentage of ground states with angular momentum L = (solid line), L = 1 (dashed line) and 
L = N (dotted line) in the vibron model with random one- and two-body interactions obtained in a mean-field 
analysis for 10 < iV < 20. 
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FIG. 7. Probability distribution P{R) of the energy ratio R obtained for A'^ = 19 and 100,000 runs for the 
spherical (solid line), deformed (dashed line) and p-boson condensate (dotted line) equilibrium configurations, 
respectively. 
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FIG. 8. As Fig. @, but for N = 20. 
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